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Abstract. For more than ten years now, many efforts have been done to 

f""*v identify and characterize nature of obstructed diffusion in model and cellular 

I lipid membranes. Amongst all the techniques developed for this purpose, FCS, 

/"" "s by means of determination of FCS diffusion laws, has been shown to be a 

fsj very efficient approach. In this paper, FCS diffusion laws are used to probe 

the behavior of a pure lipid and a lipid mixture at temperatures below and 
^~> above phase transitions, both numerically, using a full thermodynamic model, 

and experimentally. In both cases FCS diffusion laws exhibit deviation from 
^^ free diffusion and reveal the existence of domains. The variation of these 

domains mean size with temperature is in perfect correlation with the enthalpy 
C^) fluctuation. 
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1. Introduction 



Since the mosaic fluid concept of Singer and Nicolson [1 where lipids were con- 
sidered as a "sea" in which protein were embedded, the description of biological 
lipid membranes has evolved to a spatio-temporal heterogeneous mixture of its own 
components. It is now mainly admitted that biological membranes are organized 
in domains of different compositions and different size including nano, meso and 
microscopic scale organisation. Membrane heterogeneity may be of various types. 
Several lipid lamellar phases have been identified in model systems. Basically, lipid 
bilayer can exist in solid phases (s also named gel phases g), liquid disordered 
phases (Id also named fluid phases /), and liquid ordered phases (Z ) which are of- 
f*-) ten enriched in cholesterol. In complex lipid mixture, such as biological membranes, 

- * coexistence of these different phases may occur, leading to the formation of domains 

[2J. This led to the concept of raft as functional domains existing within biological 
membranes as reviewed by Simon and Ikonen [3] that highlighted their potentially 
t-H ubiquitous role in cell biology. A lot of work has been done on rafts for 10 years now 

[H 13 El [7] that led to the present consensus that rafts are heterogeneous membrane 
structures, rich in cholesterol and sphingomyelin and about 10 to 200 nm in diame- 
ter, highly dynamic in the lipid membranes of all eukaryotic cells [8] . Since this size 
is below the diffraction limit there are no direct images of their existence. More- 
over, the big discrepancy founded in size mainly depends on the technique used 
to reveal them. For example, donor quenching FRET analysis shows nanoscale 
domain formation (10 to 40 nm) in lipid bilayers with a similar composition to 
that of the outer plasma membrane, at temperature of 37 C, whereas macroscopic 
phase separation was not evident \§\. Similarly, l nanodomains can be detected 
by FRET in regions of the phase diagram in which confocal microscopy indicates 
only the presence of a single homogenous phase [TU] . Domains on length scales of 
30 to 80 nm in a region of Lo and Ld co-existence have also been detected using 
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atomic-force microscopy (AFM) [TTj . deuterium-based nuclear magnetic resonance 
(2H-NMR) and differential scanning calorimetry (TSJ, [13]. Beside their absolute 
size, since existence of domain in complex mixtures is dynamic, it seems relevant to 
use diffusion as a spatio-temporal probe of the local environment. This approach 
has been developed through Fluorescence Recovery After Photobleaching (FRAP) 
[H E3 US] or Fluorescence Correlation Spectroscopy (FCS) [HI HI US HI] and 
Single Particle Tracking (SPT)[21]. 

Each of the dynamic microscopic techniques cited above have typical advantage 
and disadvantages in respect to their timescale and statistics. For example, FCS is 
sensitive on the millisecond to second timescale, corresponding to diffusion charac- 
teristic time in lipid mixtures of a fluorescently labeled molecule through the focus 
with a waist of approximately 200nm. Recently, Wawrezinieck et al. [22] showed 
that FCS was a powerful tool for analyzing complex diffusions. They showed that 
direct fitting of the autocorrelation function was unable to discriminate between 
these complex diffusions but on the opposite, exploration at different space scale 
allowed a more detailed view of the environment structure. They developed a vari- 
able waist FCS experiment and, based on pure numerical simulation they succeeded 
in reassigning different type of FCS diffusion laws as a function of the probed envi- 
ronment such as rafts or meshwork (see Figfl]). This was furthermore confirmed by 
measuring FCS laws in living cells [53] and use recently at different space scale down 
to 50 nm [241 I20j . Even though extensive numerical simulation in given geometries 
showed a fairly accurate prediction to what occurs in very complex systems such as 
cell membranes, a level of organisation in between is still missing. Therefore, inves- 
tigations were conducted in order to further explore the ability of these quantitative 
FCS diffusion laws to distinguish phases during the transition in a well defined lipid 
mixture such as DMPC-DSPC. 
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FIGURE 1. FCS diffusion laws : FCS diffusion laws are obtained by plot- 
ting variation of the FCS diffusion time as a function of the size of the focal spot 
(or waist). Here are illustrated three typical different diffusion laws obtained by 
numerical simulation. One for free diffusion (blue line) and two for obstructed dif- 
fusion by raft (red line) or meshwork (green line). For both obstructed diffusion 
models, it can be seen that FCS diffusion laws are not linear and that their as- 
ymptotic behavior can be fitted using affinc model (For detailed explanations see 



DMPC/DSPC lipid mixtures has been studied extensively for a while. The 
thermodynamic parameters of this mixture have been deciphered by many different 
techniques including differential scanning calorimetry, neutron scattering, NMR, 
ESR, Raman spectroscopy and Fourier transform infrared spectroscopy |25j or AFM 
[25] . The structural characteristics of this mixture has also been approached by 
FRAP PS [57]. These studies have shown that DMPC-DSPC form non-ideal two- 
phases mixtures with small clusters of the minor phase in a continuum of the 
major phase. Determination of thermodynamic parameters also allowed to develop 
numerical simulations of the lipid mixture. Monte Carlo methods has been used to 
simulate the lateral distribution of each component in the pure gel or fluid phase 
of DMPC/DSPC mixture, firstly assuming one state and two components [55] and 
more recently computing a two-phases two-components triangular lattice model 
|29j . This model has been thoroughly used to predict size, shape and number of 
gel(<7)/fluid(/) clusters [30 and has recently been used to simulate FCS experiments 
on it [31]. 

The accuracy of the thermodynamic microscopic description of DMPC:DSPC 
mixtures provides the ideal model system to explore the relevancy of FCS diffusion 
laws. It is the scope of this paper to extend the previous works [281 EH EQl [31] 
into this direction both theoretically and experimentally. Here are reported FCS 
diffusion laws on both DMPC:DSPC (8:2 mohmol) mixtures and pure DMPC at 
different temperatures spanning the range of the phase transitions. These FCS laws 
shows deviations from a pure Brownian motion as expected. In order to understand 
the origin of this deviation and quantitatively characterize it, numerical simulations 
have been performed on these lipids. The results found here show that FCS diffusion 
laws are able to distinguish the existence of domains and that the mean size of the 
domains can be determined from it. A simple model is given to explain the shape 
of the FCS diffusion laws found here as a function of diffusion in a solid or a liquid 
environment. 

2. Material and Methods 

2.1. Material. The lipids l,2-distearoyl-sn-glycero-3phosphocoline (DSPC) and 
l,2-dipalmitoyl-sn-glycero-3-phospho-choline (DMPC) were purchased from Avanti 
Polar Lipids (Alabaster, AL.). They were used without further purification and 
kept at —20 ° C in chloroform:mcthanol(9:l, mol:mol) at lOOmM concentration. For 
FCS measurements lipids mixtures were labeled with 2-(4,4-difluoro-5,7-dimethyl- 
4-bora-3a,4a-diaza-s-indacene-3-pentanoyl)- l-hexadecanoyl-sn-glycero-3-phospho- 
choline(^-BODIPY FL C5-HPC)from Invitrogen (Carlsbad, CA, USA)). 

2.2. Monte Carlo simulation of the two phases two components lipid mix- 
ture. Our Monte Carlo simulations were directly adapted from the work of Sugar 
et al. [55] also described in [3T]. Briefly, the thermal fluctuations of the DMPC 
: DSPC lipid mixture was simulated using a two-state Ising type monolayer tri- 
angular lattice. Each lattice point is occupied by one acyl chain of either DMPC 
or DSPC. During the simulation, trial configuration are generated by means of six 
different elementary steps : 

- One that can be described as a phase transition step, consists in changing the 
state of a randomly selected acyl chain from gel to fluid or inversely. 

- Five that can be described as diffusion steps, consist in exchanging two neigh- 
boring molecules. In the Monte Carlo algorithm, 3 different processes can occur 
at 3 different time scales. One is diffusing in a liquid environment, the other is 
diffusing in gel environment and the last is changing its state. For lipids within a 
mixed environment (both gel and fluid) the probability of entering a diffusion step 
depends on the fraction of gel lipids chains (f c ,g) surrounding the two lipids that 
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are to enter a diffusion step. As extensiveiy described in [31] . a rate function r(/ Cj9 ) 
has been introduced in the simulation : 

AE 

where r$ is the value of r in a fully liquid environment (f c ,g — 0) and AE is the 
energy barrier needed for a diffusion step in an all-gel environment (f c . g = 1). As 
in the work of Hac et al. |31| . the phase transition step probability (r state ) was 
set to be equal to tq and AE/kT was set to 4.25 according to the following ratio : 

The Monte-Carlo simulations were implemented in C++ (Microsoft Visual C++ 
Version 6.0). They were run on a multiprocessor computer (4 AMD opterons dual- 
core 865 at 1.8 Ghz) using the Intel C++ compiler (ice vlO.O). For each type 
of simulations periodic boundary conditions were used. For the lipid mixture, 
simulations were made on 8:2 mol:mol DMPC-DSPC lipid mixtures at different 
temperatures(294, 298, 299, 300, 301, 302, 304, 306, 308, 310, 315 and 325 K). 
A 60x60 lipids lattice (60x120 lipid chains) (1 lipid=l lattice unit square) was 
simulated. For DMPC alone, a 40x40 lipid lattice (40x80 lipid chains) was used. 
10 8 Monte Carlo steps were made for each simulation except at 294 K (1.2. 10 8 ). 

2.2.1. Thermodynamic model. The thermodynamic model used here is extensively 
described in [5U] and [3T] Briefly, each lipid chains can exist in two states, gel (5) 
of nuid(/), these states being different in internal energy as well as in entropy. 
Therefore the total energy of one layer of the lipids in a given configuration C is 
given by : 

E(C) = n\E\ + n f A E f A + n B E B + n l B E l B 

_i_r,39 p99 1 n gf pgf , „// r?ff , n gg j?gg 
+n AA u AA + n AA iL AA + n AA a AA + n AB iL AB 

(o\ j-nsf rr9f j_r,fg tpfs 1 „// rpff 1 „gg pgg j_„9f pgf 1 „ff pff 

\ z ) '"-AB^AB "t" "'AB^AB "r" "'AB^AB ^ n BB rj BB + "'BB^BB "r" n BB rj BB 

where n A B and n A B are the numbers of the gel and fluid lipid chains and the E A B 
and E A B terms are the respective internal energies of the two states of species A 
and B. It is shown in [29] that the Gibbs free energy for a given state can be written 
as : 



(3) 



G = G + n A {AH A - TAS A ) + n f B (AH B - TAS B ) 

qf qf 1 qfqf\ qq qq ■ f f f f i qf qf i fq fq 

- n AA^AA + n BB U3 BB + n AB U AB + n AB u) AB + H AB U AB + n AB U AB 



Where AHa and AH B the calorimetric enthalpies, ASa = ^ A an< ^ ^^b = 
ij, — — are the respective melting entropies with T mi ^ and T m ^ B being the melting 
temperatures of the two pure components and Go the energy of an all gel lipid 
matrix. All the other parameters (w" ) are the nearest-neighbor interaction pa- 
rameters of a lipid chain species i in state a with a lipid chain species j in state 
(3. The whole set of parameter described above has been given values according to 
[3"T] since these values has been obtained on multilamellar vesicles. These values 
are summarized in table [TJ 

2.2.2. FCS diffusion laws in the simulation. In order to build FCS diffusion laws 
in the MC simulations, 20% of the total number of DMPC lipids chains were intro- 
duced as markers. During the simulation, these lipids move and generate intensity 
according to their position in the laser profile which is considered as a Gaussian 
with a given waist w (value of the radius at 1/e 2 intensity), 



±m,A 


= 297AK 




= isssj.mor 1 


Tm y B 


= 327.9£T 


^BB 


= 1474 J.mor 1 


AH A 


= 13.165fcJ.mo/- 1 


, ,99 
U AB 


= 607 J.mor 1 


AH B 


= 25.37kJ.mol- 1 


,,ff 
^AB 


= 251 J.mor 1 


AS A - 


= 44.31J.mor 1 .K~ l 


, ,gf 

^AB 


= 1548J.mor 1 


AS B - 


= 77.36 J.mor l .K- 1 


, J 9 
^AB 


= niQj.moi- 1 



TABLE 1. Monte-Carlo simulation parameters. Parameters of the Monte 

Carlo simulation of DMPC-DSPC mixture. The indices g and / correspond to gel 
and fluid state respectively. Values are the one given in 131 1 
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At each time step, the detected intensity in our simulations is computed assuming 
a Poisson distribution. The number of detected photons (n p h) for a particle at 
position (x,y) is given by a random variable following the Poisson distribution with 
parameter /3I(x,y) with (3 describing the collection efficiency of the setup [52] . 

To analyse intensity fluctuations, the normalized time autocorrelation function 
(ACF) is defined as : 

,.s ( 2 )/ x _ < n ph (t)n ph (t + r) > 

(5) 9 {T) ~ < n ph (t) >2 

where < . > represent a time average. In our simulations, the ACF is calculated 

from a generated intensity file after the whole Monte Carlo simulation. The software 

correlator used to compute the ACFs follows the architecture proposed by [33] and 

described in |32j . It has a logarithmic timescale, each channel having an individual 

sampling time and delay time. 

FCS diffusion laws are then built by correlating the intensity fluctuations ob- 
tained at different waist in the same simulation. The biggest waist being 2/3 of the 
simulated matrix. 



2.3. Image analysis. While running Monte-Carlo simulations, some images of the 
lattice were recorded at a chosen Monte Carlo step frequency, using an intensity 
code for the four different combination of the lipids states (l£P). This lead in easier 
numerical image analysis capacities for thresholding. 

Image correlation spectroscopy analyzes were conducted on thresholded binary 
images using the ICS plugin for ImageJ \http : //rsb.info.nih.gov/ij/\ devel- 
oped by Fitz Elliott, based on the work of Petersen [33] which can be found at 
http : // www.cellmigration.org /resource /imaging/ software. Using spatial im- 
age correlation spectroscopy module of the plugin generated a 2D correlogram of our 
images obtained from the MC simulation. This 2D correlogram was reduced to a 
monodimensional correlogram assuming a circular symmetry. The ID correlogram 
was arbitrary fitted by the following exponential function : 

(6) g(0 = g(0)( ai .e-^ 1 ^ + a 2 .e~^ 1 ^) + g(oo) 

where l C2 is the characteristic length of the domains. Correlograms were average 
over 50 different images at different time steps of the MC simulation at each tem- 
perature. The performance of this procedure was tested using simulated images of 
circular domains of different radius and showed error within 20 % accuracy. 

Direct morphoanalysis was also performed using the ImageJ plugin "particle 
analysis" . Analysis are conducted on binary thresholded image over a set of 50 
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different images obtained at different time of the MC simulation as for ICS. The 
"particle anaylsis" plugin is a simple binary derivative method that delimits areas 
of pixels having the same value (1 or 0) therefore giving a direct access to the 
number and area of the domains on each image. This method allows to plot the 
exact distribution in term of area and occurrence of the domains. For simplicity 
reason, it was decided to only keep the mean area of the domains defined as the 
following : 

(7) A dom = jfJ2 NaA 

A=1 

with A the value of the area (discretised in classes every 25 pixels) and N the 
number of occurence of this area within the image. 

Finally, the fractional area of gel lipids (5") defined as the following S" 



, g , ^C^CV* a« u^C ^'"»i"6 ^g ~ Sg+Sf 

within the image could also be determined at each temperature using this method, 
allowing therefore to scale the image in nm. This fraction was used to extrapolate 
the total area of the simulation and the mean area ((a) = f g -a g + fi-ai) occupied 

□ 2 ° 2 

by a lipid according the following value of lipid areas, a/ = 63 A and a g = 45 A 
for PC lipids HIES] . 

2.4. Preparation of multilamelar vesicles. 10 /d of a lOmM mixture of DMPC:DSPC 
(8:2 mobmol) labeled with j3 Bodipy Fl C5-HPC at a lipid molar ratio of 1:50 000 
were deposited on a glass coverslip previously extensively cleaned with ethanol, 
water:ethanol (7:3 vokvol) and chloroform. Lipids were dried under vacuum for 

one hour and hydrated with 700 (A of pure distilled water previously heated at 
50 ° C. Sample was then allowed to cool down for at least 30 min to the different 
measurement temperatures within the microscope chamber. 

2.5. FCS experiments. FCS experiments were performed on a home-built de- 
vice. This experimental setup has been extensively described in [24] . Briefly, it 
is based on an inverted microscope (Zeiss Axiovert 35M) with a NA=1.2 water 
immersion objective (Zeiss C-Apochromat)and a three-axis piezo-scanner (Physik 
Instrument, Germany). In order to avoid photobleaching of the fluorescent label 
during the experiment 3/iW of a CW laser at 488nm was used for fluorescence ex- 
citation. For experiments above the diffraction limit, a 30/im pinhole conjugated 
to the microscope object plane, defined the observation volume. After the pin- 
hole, the fluorescence signal is split by a 50/50 beamsplitter and focused on two 
avalanche photodiodes (Pcrkin-Elmer SPCM-AQR-13) through a 535 ± 20 nm flu- 
orescence bandpass filters (Omega Filters 535DF40). The fluorescence intensity 
fluctuations are analyzed by cross-correlating the signal of each photodiode with 
a ALV6000 hardware correlator. For FCS diffusion laws, waist of the laser in the 
object plane was modified by underfilling the microscope objective back-aperture 
using an adjustable diaphragm placed on the excitation optical path. The values of 
the observation area therefore obtained were calibrated by measuring the diffusion 
time of a 0.5 /uM Rhodamine-6G in aqueous solution at 22 ° C (D=280 /xm 2 .s _1 ). 
For each waist at each temperature of the study, at least 200 FCS measurements 
of 5s duration each were made by series of 10. 

2.6. Fitting of ACFs. For free Brownian two-dimensional diffusion in the case of a 
Gaussian molecular detection efficiency and accounting for negligible photophysics 
effects on the fluorophore (triplet state, blinking, bleaching...) the fluorescence 
autocorrelation function (ACF) is given by : 



(8) g (2 \T) = l + -— i — 

V J y \ ) -r N1 r_ 

where N denotes the average number of molecules in the observation area and 
Td the diffusion time.Td is linked to the laser beam transversal waist w and the 
molecular diffusion coefficient D by : 

2 

w 

^ T « = Id 

All the ACFs obtained experimentally using variable waist above diffraction 
limit, or on Monte-Carlo simulation were fitted using Eq. [8j 



3. Results 

3.1. FCS at variable radii and variable temperature on DMPC and DMPC: 
DSPC samples. FCS diffusion laws were acquired on pure DMPC lipid multi- 
lamellar vesicles and on vesicles made of DMPC:DSPC (8:2 mol:mol) both labeled 
with C5-Bodipy-PC at different temperatures below, within and above the phase 
transitions. Fig [2] A shows a typical correlogram obtained on pure DMPC at 297K 
with a transversal waist u>=218 /mi, Fig[2]B shows the same type of experimental 
correlogram obtained on the DMPC:DSPC 8:2 mobmol mixture at 302K with a 
transversal waist u>=210 /mi. On both part of the figure, the red line is the fit of 
the correlogram using Eq. [8] allowing determination of Td in the correlogram. The 
residual, which is plotted in the upper part of Fig [2] A and B, confirms the quality 
of the fit. Repeating experiments at different waists and different temperatures 
allows to obtain FCS diffusion laws. These are depicted in Fig[2]C and D for pure 
DMPC and DMPC:DSPC mixture respectively. 

FCS diffusion laws can be fit at different temperatures using the following equa- 
tion : 

2 

w 

(10) T d = —^+T do 

4U e ff 

with Td the experimental values of diffusion time at different waist w and Td the 
extrapolated diffusion time at w 2 = 0. These FCS diffusion laws also allow the 



determination of an effective diffusion coefficient named D e fj. Fig. 10 A and B 
at the end of this paper show the increase of this diffusion coefficient -D e // with 
temperature for both sample. 

As illustrated in Fig[3] A and B, fit of the experimental FCS diffusion laws exhibit 
a negative Td for temperatures below the phase transition (Fig [3] A, T < 297K) 
in the case of DMPC or below the second transition (gf : //) in the case of the 
DMPC:DSPC mixture (Fig[3]B, T < 310K). When the FCS laws are acquired in 
pure fluid phase (T > 310K for DMPC:DSPC mixture or T > 298K for DMPC 
alone), they exhibit T do values close to zero which shows a pure Brownian behavior 
of the diffusing molecules as expected. When plotting these Td value as a function 
of temperature, it can be seen on Fig[3]C and D that in both cases (DMPC alone or 
lipid mixture) the more the system is in the gel phase (or the more the temperature 
decrease) the more Td becomes negative. 

In order to confirm the experimental results and to understand the origin of 
a negative T do when gel domains are present in a fluid phase, it was decided to 
perform Monte Carlo simulations using a full set of thermodynamic parameters for 
DMPC and DMPC:DSPC lipid mixtures. 
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FIGURE 2. Experimental correlograms and FCS diffusion laws 

obtained on lipid mixtures. Part A and B Correlograms obtained on 
multilamellar vesicles of DMPC (Part A) and DMPC:DSPC 8:2 molimol (Part B) 
(square dots are experimental data) fitted using cq. [8] Residuals arc plotted in the 
upper part of the figure showing accuracy of the fit. f*art C and D Experimental 
FCS diffusion laws obtained for DMPC (Part C) and DMPC:DSPC 8:2 mohmol 
(Part D) at different temperatures. (Colors on line) 



3.2. FCS at variable radii on MC numerical simulations : Diffusion Laws. 

Monte Carlo simulation were made at different temperature below, within and 
above the two state transitions (gg : gf : //) of the DMPC:DSPC 8:2 (mol:mol) 
and below, within and above phase transition of DMPC alone (g : /). Fig HK 
shows a snapshot of the DMPC:DSPC 8:2 mol:mol mixture at 304K (gf state) as 
obtained from the MC simulation. Gel domains (in green) can be seen in a sea 
of fluid lipids (in red) (60x120 lipid chains). MC simulated correlograms of this 
mixture are depicted for different waist at this temperature in Fig[4t$, showing an 
increase of r^ with increasing waists as expected. As seen from experiments (see 
Fig [3]), Tci is also expected to increase with decreasing temperature at a given waist, 
this is illustrated in Fig|4p, confirming the validity of our MC simulations. At low 
temperature (T < 300K) or large waist (w > 201. u.) the correlation function at 
long time scale shows higher noise due to the limited number of events arising at 
these times. Nevertheless, the determination of Td is still valid. FCS diffusion laws 
could therefore be established here again by plotting r^ as a function of the square 
of the waist w 2 . 

As illustrated in Fig. [5] FCS diffusion laws are obtained for the different temper- 
ature of the MC simulations for both DMPC:DSPC mixture (Fig. [5|rV) and DMPC 
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FIGURE 3. Variation of experimental Td as a function of tem- 
perature. Part A and B Fit of the experimental FCS diffusion laws using eq. |lO| 
for DMPC (Part A) and DMPC:DSPC 8:2 molimol (Part B) allowing determination 
of Td a for each set of data. Part C and D r^ is plotted against temperature for 
DMPC (part C) and DMPCiDSPC 8:2 mohmol (part D). Square dots arc obtained 
value with error bar. (Colors on line) 



alone (Fig. [5)3). As found experimentally in the case of DMPC:DSPC mixture, 
above 310 K the FCS laws seems perfectly linear (Fig [5] A), accounting for a pure 
Brownian diffusion of the lipids. Below 310 K (Fig[5]A inset) all the FCS laws seems 
to be biphasic, exhibiting therefore deviations from a pure Brownian diffusion as 
expected for inhomogeneous media. Even below the first transition temperature gg 
: gl of the lipid mixture these FCS laws still shows a biphasic behavior. As shown 
in fig. [5] B the same biphasic type of FCS diffusion laws has been obtained with 
pure DMPC MC simulations below the transition temperature (296.5 K). 

MC simulated FCS diffusion laws can be fit using the equation [l0| as depicted in 
Fig. [6] A for DMPC:DSPC mixture and in Fig.|6]B for DMPC alone. 

As in the experiments, the value of an effective diffusion coefficient D e ff can be 
calculated following Eq. [lO] and is plotted as a function of temperature in supple- 
mentary material Fig[lO]C (DMPC:DSPC mixture) and D (DMPC alone). In both 
case, it is important to note that the global shape of the curve -D e // = f(T) is the 
same for simulation and experiments and is signed by a sharp increase during the 
transition from gel to fluid medium. 

It can be seen that, as in the experimental FCS diffusion laws, T& a decrease with 
decreasing temperature as plotted in Fig. |6p and D respectively for DMPC:DSPC 
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FIGURE 4. Images and Correlograms obtained by MC simula- 
tion of DMPC-DSPC 8:2 molimol mixture. Part A. Monte-Carlo 
snapshot of a 60x60 lipids lattice 8:2 mol:mol DMPC-DSPC mixture at 304K. Green 
domains correspond to gel lipids, red domains to fluid lipids. Part B. Normalized 
correlograms obtained from the Monte-Carlo simulation at 304K for different waists. 
Part C. Normalized corrclogram obtained from the Monte-Carlo simulation at a 
given waist (w— 8 l.u.) at different temperatures from gel (T < 300K) to liquid 
(T > 315K) and within the melting regime (300-ff < T < 315K). 



mixture and DMPC alone. These Td also reach values close to zero at T > 310A' 
for the DMPC:DSPC mixture or T > 298K for DMPC alone, confirming a pure 
Brownian diffusion of the system above this temperature. 

These observations confirm the validity of our experimental results on both 
DMPC and DMPC:DSPC lipid mixture. Both system exhibit deviation from a 
pure Brownian motion as soon as they are not in pure fluid state. These deviations 
can be seen using FCS diffusion laws and are signed by a negative Td ■ 



3.3. Origin of negative diffusion time at zero waist. As shown above in 
Fig[3p and D or Fig[6p and D, r extrapolated from the fit of the FCS diffusion 
laws (both experimental and simulated) using eq 10 exhibit growing negative value 



with decreasing temperature. This result is surprising since it is different from the 
one obtained by Wawrezinieck et al [22] in the case of microdomains (where t<2 is 
positive) since DMPC/DSPC mixture domains would have been expected to behave 
similarly. The results obtained here are closer to the one obtained in this previous 
study for meshworks (see FigHJor comparison). In this later case, Wawrezinieck 
et al [22 shown that the FCS diffusion law could be correctly described by the 
following equations : 
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FIGURE 5. FCS diffusion laws obtained from MC simulation 

at different temperatures. Part A FCS diffusion laws obtained from MC 
simulated DMPC-DSPC 8:2 mohmol mixture arc represented for different temper- 
ature below, within and above the melting regime of the mixture. Onset shows 
an enlargement for low waists at high temperature. These FCS diffusion laws ex- 
hibit two slopes at low temperature but starts to be linear with a null origin at 
temperature above 310K, indicating a pure Brownian diffusion process. Part B 
FCS diffusion laws obtained from MC simulations of pure DMPC arc represented 
for different temperature below and above phase transition (296. 3K). FCS diffusion 
laws exhibit here again two slopes below the phase transition and starts to be linear 
with a null origin at temperature above 298K. 



(ii) 



Where D Ti 



Td = 



AD„ 



s, 



conf 4Dn 



H-rf 



is the local diffusion coefficient, r 



7~conf ) 
domain 



if X 2 C < 2 
if X 2 C > 2 

is the characteristic diffusion 

^domain 



time inside a mesh, r con f is the confinement time defined as S con fT^ omam and S con f 
is the strength confinement that is a function of the probability to cross a barrier 
to go from one mesh to another, giving therefore an indication on the stiffness of 
the barrier. X 2 being the normalized coordinates of the waist (w) to the meshsize 



(r) and defined in the case of a square meshwork as X 2 
that for X 2 > 2 eq 11 can be rewritten as : 

„„2 



TVU : 
4r 2 



It has to be noted 



(12) 



Td 



*>conf 



4A, 



7, domain 
KT d 



(.1 ^ conf) 



In the case of the DMPC/DSPC mixture, one can simplify the system to a two 
state (gel and fluid) system with for each a typical diffusion time r? and T d exist. 
It has to be noted that at least for the MC simulation, at a given waist, T 9 d = 70r^. 
Since T 9 d is much bigger than T d , when a tracer moves into a gel environment it 
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FIGURE 6. Variation of Td obtained by simulations as a func- 
tion of temperature. Part A and B. Fit of the asymptotic part of the FCS 
diffusion law using cq |10| at different temperatures for DMPC:DSPC 8:2 mohmol 
(part A) and DMPC (part B) Part C and D. t^ is plotted against temperature 
for DMPCiDSPC 8:2 mohmol (part C) and DMPC (part D). Square dots are ob- 
tained value with error bar. A sharp increase is seen around 300K, close to the first 
meting temperature of the DMPC-DSPC 8:2 mohmol mixture (part C) and around 
296K, close to the melting temperature for DMPC. 



can appear as being strongly confined to this environment as compared to the 
same tracer into a fluid environment. Therefore, gel domains can be considered 
as barriers to free diffusion of the molecule in the liquid phase. It can be clearly 
seen from a movie generated using images step by step in the simulation (data not 
shown) that the tracer molecules move freely in the fluid environment while they 
seem to be unable to cross a gel environment at least during the duration of the 
movie (300 MC steps). Therefore, by analogy to the meshwork case described here 
above, r con / could be defined as a function of r? times the normalized surface of 
gel {Sg) that has to be crossed by the tracer molecule. 

FCS diffusion law obtained here could be described by the following simple in- 
tuitive model : 



(13) 



Td 




K(ri 



f 



S g l r a d ) 



if w 2 < aw 2 
if w 2 > ciWq 



With K and a given constants and w 2 , the value of w 2 for Td = 0. 

Eq[l3] shows that the first part of the the FCS diffusion law, below the cross over 
regime, is linked to the local diffusion constant D m i cro which is a function of the 



temperature. In the MC simulation, D Tl 



has been defined such that Di 



10Di 



By linear fit of the simulated FCS laws at small waists (w < 300i.it.), 



D m icro can be found for each temperature and it can be found that in the case of MC 
simulation D micro (325K) = 70D micro (294:K) (data not shown). If one extrapolate 
the value of D from the second part of the FCS diffusion law (300Lu. 2 < w 2 < 
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FIGURE 7. Com par ison of experimental and MC simulated Td 

to Td from eq. 14L Part A. Comparison of Td obtained by linear fit of 
the asymptotic part of the simulated FCS diffusion laws to Td obtained by cq. |14| 
as a function of temperature ( i.e. a function of normalized gel area). Part hJ. 
Comparison of r^ obtained by linear fit of the asymptotic part of the expcrimetal 
FCS diffusion laws to t^ obtained by cq. I14I as a function of temperature. Both 
arc normalized by means of K factor at 29- 



1200Lu. 2 ), then it is found that D mlcro (325K) = 160D micro (294if). This confirm 
that the shape of second part of the FCS diffusion law does contain more than only 
free diffusion. This cannot be verified for experimental laws since the first part of 
the FCS law is not accessible. 

As defined above, Td is the value of Td for w 2 = 0. From eq 13 it is found that : 



(14) 



T do =K< l {(l-70SS) 



With increasing temperatures, the relative normalized solid area (0 < S" < 1) is 
decreasing down to zero. In this study, 5™ was determined from MC simulation of 
the DMPC/DSPC mixture at different temperatures. Figure [7] shows evolution of 
Tdo with temperature for MC simulated FCS diffusion laws (Fig. [7] A) and experi- 
mental FCS diffusion law (Fig. [7] B) and the calculated values of Td as a function 
of Sg according to eq. Jl4| adjusted to the lowest temperature value of Td by K 
constant. Both figures [7] A and B shows a nice adequation of the calculated values 
with eq. [14] and the numerical or experimental obtained values confirming the hy- 
pothesis that solid domains act as barriers to free-diffusion of the tracers. Although 
the present model is rather heuristic than accurate, it suggests that the gel con- 
finement areas can be viewed as meshwork barriers for the diffusion of fluorescent 
reporters, but not as impermeable domains, or as raft-like domains. 
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3.4. Determination of domain size. Wawrezinieck et al [25] have previously 
shown that deviation from a pure linear regime in the FCS diffusion laws were a 
signature of heterogeneities in the probed environment. They also showed that the 
crossover regime in the FCS diffusion laws was a function of the relative size of both 
the probing waist and the domain size. In the case of mesh restriction to diffusion, 



Eq. 11 shows that the FCS diffusion law could be described by two linear laws 
crossing over each other close to the mesh size, therefore authorizing measurement 
of it. The FCS diffusion laws obtained from MC simulations in this study can also 
be described by two linear laws crossing over at a given waist, exhibiting a cross 
over regime which is a putative function of the domain size. 

These MC simulated FCS laws were analysed by a simple linear fit of the highest 
slope of the curves (from 3001. u. 2 < w 2 < 1200/. u. 2 ), using equation 10 since the 



first part of the curve is not accessible experimentally and therefore the cross over 
regime cannot be directly determined. Nevertheless, eq[l3]shows that the cross-over 
regime should occur at aw 2 , . Since w 2 , is the parameter that can be easily determined 
from both simulations and experiments, it seemed important to estimate in which 
extent it is related to the size of domains. 

For this purpose, snapshots of the simulations as the one depicted in Fig|4]A were 
made at different temperatures for the DMPC:DSPC lipid mixture. These images 
were analyzed by ICS and by direct space morphoanalysis as described previously. 

In the case of ICS analysis, all images between the two main phase transitions 
were analyzed using a biexponential function with two characteristics (l Cl and l C2 ) 
coherence lengths. It appeared that, in that case, l Cl varied between 0.5 and 1 (l.u.) 2 
which is a value equivalent to one chain or one lipid and indeed the smallest motif 
the image can contain. Whereas for the direct space morphology analysis, domain 
representing less than two lipids were not taken into consideration. Comparison 
of values found for each methods are made by measuring the number of lipids in 

the domains. This is directly obtained on the morphoanalysis and estimated by 

I — 2 
FCS and ICS in supposed circular domains of area irl 2 2 or 7r v 2 ° . There are no 
fundamental reasons for choosing circular domains except simplification for the 
comparison of each analysis. Indeed the snapshot of the simulations clearly shows 
that the domains are not circular nor square but have a very complex geometry, 
which has been already shown in detail by Sugar et al. |36) . 

o 2 ° 2 

Using surface area of aj = 63A and a g = 45A for PC lipids respectively in 
fluid phase or gel phase [TH |3S] allowed calculation of the mean radius of domains 
in nm. In Fig [8] are plotted the mean radius in nm of the above defined domains 
obtained by FCS diffusion laws on the MC simulation and compared to those viewed 
by ICS analysis or direct space morphoanalysis. This comparison shows that the 
mean domains radius obtained by our FCS diffusion laws are close to the one found 
both by morphoanalysis and by ICS analysis. This clearly shows that the w 2 , is an 
indirect indication of the mean radius of the domains. 

Moreover, Fig^K, that shows changes in the experimental Wq and the mean ra- 
dius (wq/2) parameters found by FCS diffusion laws as a function of temperature, 
exhibits a main change between 301 and 302 K (close to the first melting tempera- 
ture, transition gg : gf) and a smaller shoulder around 306 K (close to the second 
melting temperature, transition gf : ff). The shape of the Wq = f(T) curve is 
similar to the heat capacity profile of the DMPC:DSPC (8:2 mol:mol). This is also 
seen in FCS diffusion laws of DMPC alone (Fig^R). In that case, Wq (equivalently 
u>o/ 2) show an increased at 296K, close to the g : / phase transition temperature 
of DMPC(re/). This variation of Wq is also obtained in simulated FCS diffusion 
laws of DMPC:DSPC mixture and DMPC alone respectively (Fig|9p and Fig|9p). 
Whereas radius (wq/2) of domains viewed by FCS diffusion laws on MC simulations 
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FIGURE 8. Comparison of Wq obtained by three different meth- 
ods as a function of temperature. Three different methods (see text for 
details) were used to determine the mean size of the domains within the MC sim- 
ulation of DMPC:DSPC 8:2 mohmol mixture. These methods are compared as a 
function of temperature. Circles are the values obtained by determination of w in 
FCS diffusion laws. It has to be noted that variation of the domain size seen by 
ICS of FCS diffusion laws shows strikingly similarities with heat capacity profile 
(for comparison see Fig 8.6 p!35 of 1371 ) 



were found to be less than 5 nm, experimental FCS diffusion laws shows values of 
Wq values between 10 000 and 26 000 nm 2 for the lipid mixture, which finally give 
values oi Wq/2 higher than 50 nm and lower than 80 nm. 



4. Discussion 

Diffusion in this type of lipid mixture has been studied for a long time by var- 
ious techniques. Vaz et al. performed FRAP experiments 15J on DMPC: DSPC 
8:2 mobmol at different temperatures using NBD-DLPE as the fluorescent dye 
and found values ranging from D=0.25 /im 2 .s _1 at T=293K to D=7.5 fim 2 .s~ 1 
at T=308K. The experimental FCS diffusion laws obtained in this study exhibit 
value of diffusion coefficient between 0.12 /im 2 .s _1 at 294K to 5.6 /im 2 .s _1 at 310K. 
Indeed, the global increase of D obtained by FCS diffusion laws as a function of 
temperature is coherent with the one observed using FRAP technique. In complex 
media, FRAP measurements generally do not show total fluorescence recovery. This 
is described by fitting FRAP recovery curves with a so-called immobile fraction. 
This immobile fraction reflects the incapacity for one molecule to leave the pho- 
todestructed area during the time of the experiment and is therefore an indirect 
information on the restriction to the diffusion. Immobile fraction is ranging be- 
tween (free diffusion) and 1 (no diffusion). In their study, Vaz et al. showed that 
the immobile fraction varied from 0.52 at 293K to at 308K with a sharp transition 
in between, closely linked to the phases transitions gg : gl : 11 of the lipid mixture 
[15] . The FCS diffusion laws observed in our study exhibit non zero and negative 
r do below the second transition gl : 11. It has to be pointed out that the shape 
of the evolution of the immobile fraction as a function of the temperature of the 
lipid mixture is strikingly similar to the evolution of our r^ parameter. In a pure 
gel phase, one would expect to retrieve a free diffusion behavior of the molecules. 
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FIGURE 9. Variation of ti>g and uio/2 as a function of tem- 
perature for MC simulated DMPC-DSPC mixture or pure 

DMPC.tu is estimated from the linear fit using cq. llO of the FCS diffusion 
law as the value of w for Td — 0. Here is depicted plots of w (right scale) 
and wq/2 (left scale) as a function of temperature determined experimentally on a 
DMPC:DSPC 8:2 mohmol mixture (Part A) or pure DMPC (Part B) or numerically 
for DMPC:DSPC 8:2 mohmol (Part C) and pure DMPC (Part D). (Experimental 
point are linked by a solid curve(for -£•) or dashed curve (for w Q ) polynomial fit 
to guide the eye.) 



In this later case, the immobile fraction, as well as our r^ parameter should go 
back to the zero value. This is not the case here, neither it is the case for the 
immobile fraction in the study of Vaz et al. |15j . One simple explanation could be 
that the system is not in pure gel state at a temperature close to 294K. But the 
MC simulation using the thermodynamics of the system shows that more than 99% 
of the system is in the gel state. When this amount of gel is present in the mix- 
ture, percolation is expected to occur and to stop long range diffusion, which is the 
main diffusion component seen in our FCS experiments. Nevertheless, it has been 
suggested that even in gel state, ripple phase formation allows fast diffusion along 
line defects [38]. This would lead to fast fluctuations within the observed area and 
therefore a higher diffusion coefficient than expected. For example, Hac et al. [3T] 
also found using FCS a value of diffusion coefficient in the gel phase in the range 
of 0.05 to 0.1 fim 2 .s _1 which is again significantly higher than the one reported 
by other methods. Ripple phase formation is certainly not the only explanation of 
these high diffusion coefficients found in the gel phase since our MC simulations 
exhibit exactly the same behavior regarding the t% parameter as experiments do, 
regardless the fact that ripple phases do not exist in the simulations. Lipid state 
change could also be a source of fluorescence fluctuations within a given area and 
therefore lead to a new correlation time. The system studied here exhibits three 
different time scales : one characteristic of diffusion in a pure liquid environment 
(t^), one in a pure gel environment (t|) and finally one proportional to the rate of 
state change from liquid to solid (or inversely) for each lipid chains. It has been 
shown that this latter time has an influence on the apparent diffusion [3TJ 130] • If 
state changes occur more rapidly than the time needed for a tracer to diffuse into 
a gel obstacle, the obstacle itself will be able to fluctuate and to change his shape 
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and size. This is equivalent to diffusion in smooth obstacle where it has been shown 
that percolation has a much smaller effect than expected [39] . Finally, time scale of 
the method is also very important. FCS or FRAP experiments are conducted over 
tenths of seconds which does not authorize, at diffraction limited spots, to measure 
D lower than 1CP 4 or 10 _5 /xm 2 .s _1 . All these possibilities can explain the finding 
of a higher diffusion coefficient than the one intuitively expected and a non zero t^ 
in the pure gel phase at 294K. 

Another parameter of interest in our FCS diffusion laws is the value of the inter- 
cept of asymptotic linear extrapolation of the diffusion law with the abscissa axis. 
This parameter is empirically defined has the "zero diffusion time" characteristic 
waist (value of w 2 at t<j = 0, w 2 , — iD d ° )■ This empirical parameter has been 
analysed has a function of temperature. It exhibits a very good correlation to the 
differential scanning calorimetry curve with two maxima located at the same tem- 
perature (IK accuracy) for DMPC:DSPC 8:2 mokmol and one maximum at 296K 
for DMPC alone. Liquid or solid domains are known to be maximum in size at the 
maximal enthalpy [4"U1 |4"T] , it seemed therefore reasonable to analyse this w 2 , pa- 
rameter as being a function of the domain size. Wawrezinieck et al [2 2) have shown 
by using classical random walk simulations, that in given geometries, this parame- 
ter allows to estimate the exact size of the mesh or the domain. These numerical 
simulations results have been confirmed by analytical solving of the problem by N. 
Destainville [12] • The geometry of the system studied here is totally unknown and 
depends only on the thermodynamic properties of it. Therefore it seems very diffi- 
cult to obtain an analytical model that precisely define the size of the domains as a 
function of Wq. Nevertheless, by using ICS or morphometry analysis of the images 
obtained from the MC simulation, it is shown here that the radius of circular-like 
domains found in the images is close to the size of those estimated by means of the 
Wq parameter obtained from the FCS diffusion laws. This clearly shows that Wq 
is a good estimator of the size of the domains in the lipid mixture. In this study, 
circular-like domains are found to measure between 2.5 to 5 nm in radius in the case 
of MC simulations while they exhibit size in between 50 to 80 nm in experiments. 
AFM studies on a 50:50 mokmol DMPC:DSPC supported bilayers lipid mixture 
has shown to exhibit domain size between 50 (circular) and 150x300 nm (rectan- 
gular) in the g : / phase coexistence [43]. These values are coherent with the one 
observed here experimentally for our 80:20 mokmol DMPC:DSPC lipid mixtures, 
diss et al. [H] also showed existence of nanoscale domains in DMPC-DSPC mix- 
ture. AFM imaging allowed them to measure gel domains of 10 nm, while neutron 
diffraction permitted identification of gel domains exhibiting values of 7 nm. These 
size are comparable to the one observed here in our MC simulations whether by 
ICS analysis or by estimation of the w 2 parameter in the simulated FCS diffusion 
laws. 

As stated above, this two-phases two-components DMPC:DSPC lipid mixture 
has been studied for a long time by many different techniques both on dynamic and 
structural approaches. Bagatolli et al. [45] even showed existence of micrometer 
size domains within GUVs made of DMPC:DSPC 1:1 mokmol. 

Indeed, all these observations show that such systems exhibit different space 
and time scale confinement, rising the question of the anomality of the diffusion. 
Anomalous diffusion has been extensively used to describe and study many differ- 
ent type of dynamic behavior of biological molecules in their complex media (for 
review see [46]). Basically, it can be seen as a general case for diffusion. Relation 
between time and space could be generalized as < r 2 >cx t a , with < a < 2. If 
a = 1, the diffusion is normal and described by pure Brownian motion whereas 
if a < 1 the system is sub-diffusive indicating that the molecular motions are 
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restricted. DMPC:DSPC two-phases two-components lipid mixture has been de- 
scribed by means of anomalous diffusion both on MC simulation [30] and experi- 
mentally [31 j - In this case, a was found to be systematically less than one, except in 
the / : / region and in the far g : g region (very low temperature) , where the molec- 
ular motion becomes Brownian again. FCS diffusion laws obtained on both MC 
simulations and experiments have also been fitted using anomalous diffusion. In this 
study, a has been found to vary between 0.64 and 1 in both cases (MC simulations 
and experiments) exhibiting two minima at the phase transitions g : g — > g : f and 
g :/—>/: / (data not shown). Anomalous diffusion can therefore clearly describe 
the molecular motion in this two-phases two-components lipid mixtures or during 
phase transition in a pure lipid, nevertheless, even if it gives informations on the 
heterogeneity of the systems, it fails to quantitatively interpreted the structure of 
this heterogeneous medium. It has been a deliberate choice all along this work to 
compare results at extremely low spatial scale (MC simulation) with experimental 
scale, this later being 3 orders of magnitude larger. This is motivated by anoma- 
lous properties of diffusion in the lipid mixture that exhibit similar space and time 
behavior at any scale. This question of time to space invariance can also be viewed 
a posteriori by the consistency of the results found here (wq versus T for example) 
although simulations and experiments are compared at very different scales. 

5. Conclusion 

Here is presented a study of the two-phases two-components DMPC:DSPC 8:2 
mol:mol lipid mixture and pure DMPC lipid system by means of FCS diffusion laws 
both experimentally and using MC simulation with a complete thermodynamic de- 
scription. It is clearly shown that these FCS diffusion laws, allows quantitative 
characterization of the system in terms of diffusion, phase transition and mean size 
of the gel or fluid domains presents in the lipid mixture. Domain size and transient 
confinement times has already been shown to be predictable by the use of FCS 
diffusion laws on defined geometries and pure random walk model |22) . Here it is 
shown that this can be generalized to a more complex model where the geometry is 
unknown and the molecular motion are only driven by the thermodynamic param- 
eters of the system itself. Finally, the parameter obtained from the interpretation 
of the FCS diffusion laws issued of the numerical simulation are shown to have the 
same behavior than the one obtained on experimental system. This confirms the 
advantages of using FCS diffusion laws on experimental systems to described their 
temporal and spatial structure. 
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FIGURE 10. Effective diffusion coefficient obtained by FCS dif- 
fusion laws. Variation of the diffusio n co efficient (D e ff) obtained by measuring 
the slope of the different fits using cq. |10| of the FCS diffusion laws as a function 
of temperature. This clearly show, as expected, an increase in D with temperature 
from a lower plateau (gel phase) to a higher one (liquid phase). Part A and B arc 
values experimentally obtained whereas part C and D are numerically obtained for 
DMPCrDSPC 8:2 mohmol and pure DMPC respectively. 
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